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I. INTRODUCTION 

In fusion experiments, it is a common method to heat the plasma by means of microwaves 
resonating with the electron cyclotron frequency Uce or its harmonics. The underlying mech- 
anism is well understood^, it allows to heat the electrons and drive significant currents. If, 
however, the plasma density exceeds the cutoff density for the frequency of the injected mi- 
crowave, it is reflected before reaching the resonance and heating the plasma. This problem 
can be overcome by heating at higher harmonics of the cyclotron frequency, which requires 
electron temperatures of several keV to be efficient^. Another approach is the utilization of 
electron Bernstein waves (EBWs). These are electrostatic waves that cannot propagate in 
vacuum and therefore need to be coupled to electromagnetic waves via mode conversion. No 
density cutoff exists for EBWs, and they are very well absorbed at Uce and its harmonics, 
even for low temperatures. Details about the physics and applications of EBWs can be 
found in a recent review article^. 

EBW heating has been successfully demonstrated in stellaratora^ and tokamaks^. Es- 
pecially in spherical tokamaks, there is a need for starting and sustaining a plasma with 
reduced induction current since their geometry leaves only little space for a central trans- 
former coil. Therefore, other current generation mechanisms have to be applied during the 
current ramp-up phase to save magnetic flux from the small solenoid. EBW current drive 
is one of such mechanisms, as demonstrated in TST-2^ and MAST-. EBW heating is also 
useful during the ramp-up phase, as it increases the electron temperature and decreases the 
collisionality, thus resulting in more effective induction. Helicity injection^*^ and fast wave 
heating and current driveiS are also expected to benefit from electron heating by EBWs in 
Pegasus. Synergistic combinations of these start-up techniques have the potential to create 
a suitable plasma target to hand over to Ohmic heating and neutral beam injection in larger 
spherical tokamaks including a possible future component test facility^. 

In NSTX, the effect of collisions on the propagation of EBW has been studied experimen- 
tally^ and numerically^'^. For the Pegasus Toroidal Experiment, numerical studies of 
EBW propagation from the conversion layer and damping at the outermost Doppler-shifted 
electron cyclotron resonance (ECR) have been performed in the past, mostly by means of 
the GENRAY ray tracing code and the CQL3D Fokker-Planck code^^. This paper investi- 
gates the coupling of externally launched electromagnetic waves with EBWs, using the 2D 
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full-wave code IPF-FDMOi^. The full-wave approach is indispensable to properly model 
the propagation of the O- and X-waves and their conversion into EBW because the vacuum 
wavelength of the injected microwaves, Aq ~ 12.2 cm, is comparable with the size of the 
plasma (average plasma minor radius a ~ 40 cm). 

The efficiency of the coupling is studied as function of the toroidal and poloidal injection 
angle. The role of the shape of the plasma density profile in front of the antenna is investi- 
gated and the possibility of heating at higher harmonic of the ECR frequency is explored. 
Thus, essential information for the design of an EBW heating system can be gained. 

This paper is organized as follows: in section [Tll the Pegasus Toroidal Experiment 
is introduced and the equilibrium used in the modeling is described. The full-wave code 
IPF-FDMC is described in section IIIII Plasma and wave parameter scans in ID and 2D 
simulations are presented and discussed in sections [IV] and |Vl respectively. We scan the 
density profile, the vertical position of the plasma, the injection angles, the polarization and 
the frequency of the injected microwave. Section IVTl summarizes and concludes this paper. 

II. THE PEGASUS TOROIDAL EXPERIMENT 
A. General properties 

The Pegasus Toroidal Experiments^ is a low aspect ratio facility with values of A = 
R/a = 1.15 — 1.3, where R and a are the major radius and the average plasma minor radius, 
respectively. It was designed to explore the operational regime for spherical tokamaks with 
high toroidal beta plasmas in the limit of A — 1^". EBW heating can assist during the start- 
up phase of Pegasus as it has been done in MAST—. The plasma core can be heated-up and 
EBW heating can provide a tool to control the current profile. Furthermore, as discussed in 
the introduction, synergies are possible between EBW heating and helicity injection^i^ and 
between EBW heating and fast wave heatingi^. 

EBW emission, heating and current drive experiments are under consideration for Pega- 
sus. An antenna would be installed for this purpose in one of the 12 equatorial ports. The 
diameter of the port (40 cm) puts a constraint on the maximum diameter of the injected or 
collected microwave beam. Figure [1] shows various Pegasus' flux tubes and the orientation 
of the simulation plane. 
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FIG. 1. Pegasus flux surfaces, including a field line on the q = 5 flux surfaces, with q being 
the safety factor. The simulation plane is indicated on the right, the antenna is located in the 
equatorial plane at i? = 1 m. 

Typical magnetic field strengths are on the order of 0.1 T in the core, making 2.45 GHz a 
good candidate for central EBW heating at the fundamental harmonic. Higher commercial 
and industrial standards, such as 3.6 and 5.5 GHz, are also of interest, for 2^*^ harmonic 
EBW heating in the present configuration, or for 1^^ harmonic heating in a possible Pegasus 
upgrade to higher magnetic fields. Typical densities in the plasma center are on the order 
of He = 10^^ m~^ and electron temperatures up to Te ~ 300 eV are achieved. 

B. Equilibrium used in the modeling 

The KFIT equilibrium code^*^ was used to model high performance plasma targets (of 
plasma current Ip = 150 — 300 kA) for EBW experiments. Contour plots of the three 
components of the magnetic field in the poloidal plane as obtained from KFIT are shown in 
Fig. El Microwaves will be injected from the position R = 1 m and Z = m. 

A reasonable plasma pressure profile inside the last close flux surface (LCFS), which, 
for this equilibrium, is located at i? = 0.72 m for Z = m, serves as input constraint for 
the KFIT code. With the corresponding density profile and the values of the magnetic field 
strengths, the characteristic cutoff and resonance frequencies can be calculated. They are 
shown in Fig. |3] as function of the radial coordinate R. The value of the electron cyclotron 
frequency fee is below the values of the cutoff frequencies in the complete cross section. 
Hence, the cyclotron frequency is shielded by cutoffs, and the entire Pegasus plasma from 
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FIG. 2. Pegasus magnetic field adopted for the full-wave modeling from the KFIT equilibrium 
code. Shown are, respectively, the contours of the radial, vertical and toroidal component in a 
poloidal cross section. The values of the magnetic field strength are given in mT. 
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FIG. 3. Cutoffs and resonances inside the LCFS as function of the radial coordinate for Z = m, 
different line styles represent different cutoffs and resonances, as labeled in the plot. 

the LCFS inwards, can be referred as overdense. It is also visible that the cutoff frequency 
at the LCFS is still well above 2.45 GHz. Thus, the 0-mode cutoff layer (and potential 0-X 
mode conversion layer) for 2.45 GHz lies well outside the LCFS. In the mode conversion 
region, in the absence of density data from KFIT, assumptions about the shape of the 
density profile are necessary. These assumptions are very delicate because the shape of the 
profile, and in particular the density gradient length in the mode conversion region, play a 
crucial role in the mode conversion efficiency. Due to this sensitivity issue, different profiles 
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FIG. 4. Radial density profile at Z = m as used in the full-wave modeling: (a) assumed core 
profile and (b) edge profiles assumed on the basis of Langmuir probe measurements with the cutoff 
densities for different microwave frequencies marked by horizontal lines. 

were used in modeling the mode conversion region. Figure UK shows the assumed radial 
density profile at Z = m. This is the internal density profile valid inside the LCFS. It 
was extrapolated outside the LCFS, up to the antenna, located at -R = 1 m. Three different 
profiles, which are based on preliminary Langmuir probe measurements recently performed 
in the scrape-off layer, were used for the simulations presented in this paper and are shown 
in Fig. IDd: 

■ P 



ne,o exp <^ - — 



ne,3{p) = ne,i{p) - hi tanh <^ h 



P - P2 

W2 
P - P3 
W2 



(1) 
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where p is the normalized radius, Uefi = 2 ■ 10 m""^, wi = 0.7, a = 3.1, p2 = 1-37, ps = 1.48, 
6i = 4 • 10^^ m~^ and 62 = 16. In the experiment, it is speculated that such different shapes 
can be actively realized by different positions of a limiter placed closed to the antenna (i. e. 
outside of the LCFS). 



III. THE FULL- WAVE CODE IPF-FDMC 

From Fig. Hb, it can be deduced that, in Pegasus, the density gradient length scale 
Ln, normalized to the vacuum wavelength, can become as small as k^Ln ~ 2 in the mode 
conversion region. For these steep profiles, with respect to the wavelength of the injected 
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microwave, the geometric optics assumptions are not valid. It is therefore not possible 
to estimate the 0-X-B mode conversion efficiency by analytic formulas^""— based on the 
validity of these assumptions. Hence, a full-wave description is necessary to model the mode 
conversion process. The code IPF-FDMC provides such a modeling. 

In this section, the conversion process that results in the generation of an EBW is briefly 
described followed by a description of the full-wave code IPF-FDMC. Results of modeling 
the mode conversion process in Pegasus are described in sections lIVI and IVl 

A. O— X-B mode conversion scheme 

The 0-X-B mode conversion scheme was conceived in 1973 by Preinhaelter and Kopeckyi^ 
as a method to heat overdense plasmas. In this scheme, an 0-mode needs to be injected 
from the low field side at an optimum angle with respect to the background magnetic field. 
This 0-mode is converted into an X-mode in the vicinity of the 0-mode cutoff layer, which 
then propagates outwards until it reaches the upper-hybrid resonance (UHR) layer. There, 
it converts into an EBW, propagating backwards, which can then be absorbed at the ECR 
and its harmonics. Non-optimum injection angles lead to less efficient conversion, due to 
partial reflection from the cutoff layer or, for strongly non-optimum angles, refraction in the 
non-uniform plasma. 

B. The full- wave code IPF-FDMC 

The full-wave code IPF-FDMC is a finite-difference time-domain (FDTD) code, solving 
Maxwell's equations, coupled with an equation for the current density J, which is obtained 
from the fiuid equation of motion of the electrons, on a Cartesian grid. Details on the FDTD 
method in general can be found in Ref.'^. The system of equations to be solved for each 
time step reads: 





-V X E 



(6) 



(4) 



(5) 
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where E and B are respectively the electric and magnetic field of the wave, Upe is the elec- 
tron plasma frequency, Uce the electron cyclotron frequency and u a collisional frequency. 
The z/J term is responsible for collisional damping of the wave. While this damping mech- 
anism is real and non- negligible, here it is artificially enhanced to prematurely cause the 
complete damping of the wave before its wavelength becomes too short (when the X-wave 
is approaching the UHR) for the adopted numerical grid to resolve it. Note that this en- 
hancement does not change the conversion efficiency^^. A more realistic picture is that the 
slow X-wave experience some collisional damping at the UHR and converts into an EBW 
that propagates towards the plasma core and is cyclotron-damped in the vicinity of the out- 
ermost ECR. This treatment, however, would require a finer numerical grid, finite Larmor 
radius corrections and a model for the cyclotron-damping of EBWs suitable for full-wave 
calculations or, alternatively, the coupling of the 0-X full-wave solution with a ray tracing 
code for the propagation and damping of EBWs^i. Such improvements are left for future 
work. 

Here, B, E and J are three dimensional vector fields and all three components in the 
R, Z and direction are calculated and advanced in time, although they are treated as 
functions of R and Z only. In other words, the plasma inhomogeneity is 2-dimensional 
and the problem is invariant under translation in the transverse direction. The code was 
successfully used to model the 0-X-B mode conversion process in the TJ-II stellarator^- 
and the RFX-mod reversed field pincb^. Furthermore, microwave heating (although not 
by EBW) of the TJ-K stellarator— has been modeled with IPF-FDMC and a predecessor 
of the code has been applied to the WEGA stellarator More details about the code are 
given in Ref.— . 

Equations (jll)-(l6]) include only cold plasma effects, which are sufficient to model the 
0-X mode conversion. In the vicinity of the UHR, the X-mode becomes more and more 
electrostatic, its wavelength becomes shorter and the cold plasma formulation breaks down. 
To resolve this singularity, it is in principle possible^^ to include first order finite Larmor 
radius corrections^^ in the code to account for the X-B conversion. However, here, electron 
collisions are used to damp the X-mode around the UHR, since the inclusion of the X- 
B conversion would significantly increase the computational time. The reason is that the 
wavelength of the EBW is comparable with the electron Larmor radius, i. e. much smaller 
than the vacuum wavelength, thus requiring a much finer numerical grid. Therefore, we 
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restrict ourselves to the 0-X conversion in this paper. Note that some effects that deteriorate 
the overall conversion efficiency are included in the present model, but some others are not. 
A deteriorating mechanism is the coupling between internal slow X-waves (excited via the 
0-X conversion) into external fast X-waves, leaving the plasma^i. Slow and fast X-waves 
are normally separated by an evanescent layer which, however, can be very thin in the low 
magnetic field edge of a spherical tokamak. This effect is automatically taken into account by 
the full- wave code. An effect, not taken into account by the code, is the parametric decay 
in waves of different, generally undesired frequencies^. The effect of density fluctuations 
on the 0-X conversion, which are expected to locally and temporarily modify the optimal 
direction for efficient mode conversion, is also not taken into account, but will be the subject 
of a future study. On average, this effect is a reduction of the conversion efficiency. Hence, 
the 0-X conversion efficiency deduced from the present simulations has to be considered as 
an upper limit for the actual overall 0-X-B conversion efficiency. Note also, however, that 
the X-B conversion does not introduce any additional angular dependence. Therefore, the 
optimal direction and angular tolerance of the actual 0-X-B process coincide with the 0-X 
optimal direction and angular tolerance obtained from the present simulations. 

The vacuum wavelength of the injected microwave is set to be 256 grid points on the 
numerical grid and the normalized collision frequency lies in the range 10~^ < v/ujq < 10~^ 
(see Ref.— for details). 

IV. ID SIMULATIONS 

Calculations on a ID grid, which require significantly less computational resources than 
calculations on a 2D grid, were performed first. The plasma density is taken to be a function 
of the radial coordinate only, and the magnetic field is taken to be perpendicular to the 
radial coordinate with a constant value of ujce/^Q = 0.6. The injected microwaves correspond 
to plane waves in this case. Although the experimental situation is not described properly 
by these simplifications, the results from these calculations can serve to set constraints on 
parameters for the 2D calculations. 

For small values of koLn, Mj0lhus' formula for the conversion efficiency as a function of the 
injection ang is no longer valid. This is nicely illustrated in Fig. |5l where simulations for 
three different values of koLn are compared with the corresponding solution from Mj0lhus' 
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FIG. 5. Conversion efficiency as a function of tlie parallel component of the normalized refractive 
index, where parallel refers to the direction of the magnetic field. The different symbols corre- 
spond to results obtained from full-wave simulations for different values of the density length scale 
normalized to the vacuum wavelength, koLn, as labeled in the plot. The solid lines indicate the 
corresponding solution from the Mj0lhus equation^!. 

equation. For these ID simulations, the density profiles were taken to be of parabolic shape. 
In the case considered here, the variation between the analytical solution in the WKB limit 
and the full- wave solution is not very strong for koLn > 5. However, at small density gradient 
length, for example such that k^Ln = 2, the WKB limit and Mj0lhus' formula are no longer 
valid. This results in the discrepancy between the values obtained from the simulations and 
Mj0lhus' formula that can be seen in Fig. [51 The general trend that for steeper profiles the 
conversion efficiency becomes less sensitive to an angular mismatch can also be clearly seen. 

If the density profile becomes too steep, the polarization of the injected microwave needs 
to be re- adjusted to obtain maximum conversion efficiency^^. This effect is illustrated in 
Fig. [6], where the optimum polarization Er/ E±_ (with Er the wave electric field along the 
calculation grid and Ex_ the field perpendicular to R and to the direction of the magnetic 
field) is plotted as function of k^Ln- The values were obtained from a series of simulations 
in which the polarization was scanned in order to find its optimum value. As one can see, 
for values of k^Ln < 5, the polarization differs from its asymptotic value, which can be 
calculated analytically^ for this configuration to Efi/E± ^ 0.75. 

To check the relevance of the aforementioned polarization adjustment, two simulations 
were performed for steep profiles with koLn = 2, one with the adjusted polarization and the 
other with the non- adjusted polarization. Their results are shown in Fig. [71 One can clearly 
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FIG. 6. Polarization of injected microwave beam, Er/E± with Eji the wave electric fields along the 
calculation grid and E± the field perpendicular to R and to the direction of the magnetic field, for 
optimum conversion efficiency, obtained from full-wave simulations, as function of the normalized 
density gradient length. 
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FIG. 7. Conversion efficiency as function of the parallel component of the normalized refractive 
index obtained from simulations with koLn = 2 with optimized and non-optimized polarization of 
the injected beam. 



see that optimizing the polarization has only a fairly small effect on the conversion efficiency 
and on the optimal injection angle, or, in other words, the conversion efficiency is relatively 
stable against polarization mismatch. 
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V. 2D SIMULATIONS 



From the ID model discussed in the previous section, it could be deduced that for steep 
density profiles, the angular window to achieve maximum mode conversion efficiency is fairly 
large. That result is confirmed here in 2D, taking into account the geometry of the plasma 
and the finite size of the beam. 

The flux surfaces in Pegasus are significantly curved, as can be seen in Fig. [1] Their 
curvature radius in the region of interest is on the order of 50 cm and non-negligible on the 
transverse length scale of a realistic 2.45 GHz beam. A previous work^^ showed that in the 
presence of curved flux surfaces it is important for wavefronts to match their curvature to 
the curvature of the conversion layer. Hence, a Gaussian microwave beam emitted from the 
antenna with the beam waist located inside the plasma is considered in the simulations (for 
details on Gaussian beam properties, see e. g. Ref.~). Following a previous, internal design 
study, the emitting microwave antenna is located at i? = 1 m with the beam waist located 
at -R = 0.7725 m and a beam radius (at the waist) of wq = 0.67Ao. 

In the 2D simulations, the poloidal and toroidal injection angle can be scanned in order to 
find the maximum conversion efficiency. Figure M shows three snapshots at different times 
of the absolute value of the wave electric field for a poloidal and toroidal injection angle 
of 10° and 0°, respectively. Density profile #1 (see Fig. |Dd and Eq. ([1])) was used in this 
simulation. The positions of the characteristic frequency layers for 2.45 GHz are given in 
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the plot. In the first snapshot, taken after T = 2 oscillation periods, the wave has not yet 
reached the conversion layer. Its focused phase fronts can be clearly seen. At T = 3.5, the 
wave has reached the conversion layer and one can see a complicated interference pattern 
building up between the incoming and the reflected (i. e. not converted) wave. After T = 20 
oscillation periods, a steady state situation has been reached and a conversion efficiency of 
33% is found. The generated X-mode, which is visible by its enhanced wave electric field 
and the small scale structure, is damped at the UHR, since no X-B conversion is included 
in this simulation (see Sec. IIIIB|) . 

In order to check for the infiuence of the shape of the density profile on the conversion 
efficiency, for each of the three profiles shown in Fig. |Dd, more than 500 runs have been 
performed, each with a different combination of poloidal and toroidal injection angle. The 
result of these calculations is shown in Fig. |9l where the conversion efficiency is plotted as 
a function of the toroidal and poloidal injection angle. The maximum conversion efficiency 
is found at different injection angles for different density profiles. This is especially true 
for profile #1, as opposed to the two other profiles. The general shape of the angular 
window differs also when comparing the results from profile #1 with the two others: the 
window becomes broader with the density profile moving closer to the antenna. The highest 
conversion efficiencies of 75% are also achieved for profile #3, in which the distance between 
the mode conversion layer and the antenna is the shortest. Density profile 7^2 and 7^3 have 
shorter gradient lengths in the mode conversion region, and, thus, smaller values of k^Ln, 
which results in the observed broader angular window (c. f. Fig. |5]). 

These results clearly illustrates the importance of the knowledge of the actual density 
profile in front of the antenna for efficient EBW coupling. It was also shown how smaller 
values of koLn result in larger angular windows for efficient mode conversion (the same result 
was found in the ID simulations). Hence, steeper profiles in front of the antenna seem to 
be preferable in this case, although even steeper profiles would lead to a deterioration of the 
conversion efficiency, as described in Sec. IIIIBI 

A. Stability against vertical displacement of the plasma 

Highly elongated plasmas like Pegasus can be subject to vertical displacements during 
the discharge. It is therefore of interest, how an EBW heating system couples with a plasma 
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FIG. 9. 0-X conversion efficiency as a function of the poloidal and toroidal injection angle for (a) 
density profile #1 (b) #2 and (c) #3 (see profiles in Fig. 

that moves vertically. Thus, simulations have been performed, where the antenna was shifted 
10 cm upwards and 10 cm downwards relative to the plasma, corresponding to shifting the 
plasma 10 cm downwards and upwards, respectively. 

Figure [TU] shows the 0-X conversion efficiency for the three density profiles with a dis- 
placement of +10 cm, without and with —10 cm displacement. Note that such vertical 
displacements are extreme and rarely observed. Even under these extreme circumstances, 
the optimal launch angle in the vertical direction only varies by approximately 5° for the 
density profile #1, i. e. by much less than the angular window width. For the two other 
profiles, the variation is slightly stronger, but still below 10°. Hence, the plasma is relatively 
stable against vertical displacement and a reasonable amount of the injected microwave 
power would still be converted into an EBW when a small vertical displacement occurs. 

B. Mode conversion efficiency at higher microwave frequencies: 3.6 and 
5.5 GHz 

One feature of EBWs is that they are very well absorbed at harmonics of the ECR 
frequency, in contrast to conventional ECR heating, which requires electron temperatures of 
several keV to be efficient for higher harmonic heating^. For this reason, simulations with a 
microwave frequency of 3.6 GHz have been performed which corresponds to 2^*^ harmonic 
heating. In Fig. [TTk-c. the obtained conversion efficiencies for the density profiles 7^1-#3 are 
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FIG. 10. 0-X conversion efficiency as a function of the poloidal and toroidal injection angle for 
(a)-(c) density profile #1 with the plasma shifted 10 cm upwards, not shifted and shifted 10 cm 
downwards, respectively, (d)-(f) the same for density profile #2 and (g)-(i) for density profile #3. 

shown, respectively. When comparing them with Fig. |9K-c, where the conversion efficiencies 
for 2.45 GHz are shown, one can see that slightly higher values (75-80%) are achieved at the 
optimum injection angles, which, for their part, have barely changed. The angular window 
for high efficiencies has become smaller, which is due to the fact that with the increased 
value of fco ~ 75.5 m~^ (/cq ~ 51.3 m^^ for 2.45 GHz), the value of fco-^n at the conversion 
layer has also increased, which results in a stronger sensitivity to the injection angle, as 
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FIG. 11. 0-X conversion efficiency as function of poloidal and toroidal injection angle for (a)-(c) a 
microwave frequency of 3.6 GHz and the density profiles #l--#3 and (d)-(f ) a microwave frequency 
of 5.5 GHz and the same profiles. 

discussed in Sec. HVl 

A possible upgrade of PEGASUS to higher magnetic fields would require higher microwave 
frequencies, due to the increased ECR frequency. To check for the potential of such an up- 
grade on EBW heating, additional simulations were performed for 5.5 GHz. The conversion 
efficiencies obtained are shown in Fig. [TTH-f. One can see that the angular window has 
become slightly smaller due to the increased value of k^Ln- The difference between the 
shape of the contours of the conversion efficiency for the three different density profiles is 
smaller as compared with the lower microwave frequencies. This is due to the reduced rela- 
tive difference of k^Ln between the different profiles. The maximum efficiencies achieved for 
5.5 GHz are similar to the case for 3.6 GHz. Thus, this scenario seems also to be suitable 
for EBW heating. 
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VI. SUMMARY 

The 0-X mode conversion process has been modeled with the full- wave code IPF-FDMC 
in the Pegasus Toroidal Experiment. A frequency of 2.45 GHz has been chosen, since 
experiments are under consideration at this frequency. Different density profiles, based on 
preliminarily Langmuir probe measurements, were included in the simulations. All of these 
profiles are rather steep, with values of the normalized density gradient length of koLn ~ 2 
in the mode conversion layer. Maximum conversion efficiencies on the order of 70-75% are 
found. The angular window for conversion efficiencies above 50% is fairly large for all density 
profiles. However, they are of different shape, which illustrates that the knowledge of the 
shape of the actual density profile in the mode conversion region is important, even for such 
small values of koLn, if an EBW heating and current drive system with a power level in the 
MW regime is considered. 

An extreme vertical displacement of the plasma by ±10 cm, corresponding to an extreme 
instability, still results in high conversion efficiencies on the order of 65%, if the injection 
angles are not adjusted. Hence, the 0-X conversion is relatively stable against such dis- 
placements, thanks to the large angular windows mentioned above. 

For 2^*^ harmonic heating with 3.6 GHz, higher conversion efficiencies of 75-80% are 
obtained, but the width of the angular window is slightly smaller, due to the increased 
value of koLn at the mode conversion layer. A potential upgrade of PEGASUS to higher 
magnetic field strengths would require higher microwave frequencies, such as 5.5 GHz, for 
which simulations yield similar conversion efficiencies as for the case with 3.6 GHz, but with 
yet smaller angular windows. 

To conclude, the full-wave simulations showed that high 0-X conversion efficiencies of 
up to 80% can be achieved at Pegasus. These estimates represents an upper limit for 
the overall 0-X-B conversion efficiency, the reason being that some degradation effects like 
the excitation of parametric instabilities during the X-B conversion are not included in the 
simulations yet. It should also be pointed out, however, that the X-B conversion does not 
introduce any further dependence on the launch angles. Therefore, the optimum angles for 
0-X conversion calculated in this work are also the optimum angles for EBW heating and 
current drive at Pegasus by means of the 0-X-B conversion. 
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